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Abstract. We theoretically investigate the magnetic-field-angle dependence of the flux-flow 
resistivity pi in unconventional superconductors. Two contributions to pi are considered: one is 
the quasiparticle (QP) relaxation time r(ky) and the other is ljoCsf), which is a counterpart to 
the interlevel spacing of the QP bound states in the quasiclassical approach. Here, fcp denotes 
the position on a Fermi surface. Numerical calculations are conducted for a line-node s-wave 
£_) ■ and a d-wave pair potential with the same anisotropy of their amplitudes, but with a sign 

change only for a d-wave one. We show that the field-angle dependence of pt differs prominently 
between s-wave and d-wave pairs, reflecting the phase of the pair potentials. We also discuss 
the case where r is constant and compare it with the more general case where r depends on k^ . 

G\ 
<N 
_t ■ 1. Introduction 

©The Cooper pairing mechanism is reflected by the symmetry of the superconducting pair 

potential. Therefore, the elucidation of the pair potential symmetry is of great importance 

for obtaining the clue to the Cooper pairing mechanism in unconventional superconductors [lj. 

There are two factors of the pair potential as a complex number: one is the amplitude and the 

other is the phase. The magnetic-field-angle resolved thermal conductivity and specific heat 

measurements are powerful techniques which can detect the anisotropy of the pair potential 

amplitude [2j [3]. However, they cannot detect the phase of the pair potential. It is crucial 

to probe its phase in order to discriminate unconventional sign reversed pair potential from 

conventional sign conserved one. To this end, we propose that the field-angle dependence of the 

flux-flow resistivity pf [I] can be a phase-sensitive probe. 

The flux-flow resistivity is pt(T) oc r(e = koT) in moderately clean systems |5]. £ is the 

energy of the quasiparticle (QP) bound state inside a vortex core and T is the temperature. In 

a previous paper [6j, we took into account only the QP scattering rate r(kp) as a contribution 

to pf. Here, the Fermi wave number kp denotes the position on a Fermi surface (FS). Actually, 

there is another contribution to pf, which is an energy scale wo(^f) related to the vortex bound 
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state spectrum [JJ. In this study, we investigate the field-angle dependence of pf taking into 
account the contribution of a;o(fcF) in addition to r(k F ). 

2. Formulation 

@We consider a single vortex at low magnetic fields. We focus on the QP scattering due to 
the non-magnetic impurities distributed randomly inside a vortex core. The QP scattering 
rate r is obtained by calculating the imaginary part of the impurity self energy ImZ 1 , which 
corresponds to the energy level width 5E. The QP scattering rate r(k F ) is related to 5E as 
r(k F ) ~ ImI7 ~ SE ~ l/r(fcp). In this paper, we set fi = 1. On the basis of the quasiclassical 
approximation [8] and the Kramer-Pesch approximation [9], the quasiclassical Green's functions 
and the impurity self energy in the vicinity of a vortex core are obtained analytically |10j . 
From the pole of the regular Green's function, the QP energy spectrum E(y, k F ) is obtained as 
E(y,k F ) = 2y\d(kF)\ 2 AQ/vF±. Here, y is the impact parameter, d(kp) is the anisotropy factor 
of the pair potential, Ao is the maximum amplitude of the pair potential in the bulk, and v F ± 
is the FS average of |f f_l(^f) I defined in ref. [6j. In the quasiclassical approximation, the QP 
spectrum is continuous with respect to the impact parameter y [llj. The impact parameter 
y(> 0) is related to the angular momentum I of the QPs running around a vortex core such 
that — / = y|fcF±|- Here the negative angular momentum is due to the direction of the circular 
motion of the QPs around a vortex core. The spectrum of the low-energy vortex bound states 
is expressed [121 03] as E{y,k F ) = —uJo(k F )l = coo(k F )y\k F ±\ in terms of a counterpart to the 
interlevel spacing of the QP bound states uo(kp). Comparing the above two expressions of 
E(y,kp), coo(k F ) is represented as 

2|d(fc F )| 2 Ag 
\k F± \v F ± 

Then, pf is given as [7J 

Pi (T, a M ) oc - — — = 2 ^ (2) 

{ujo(k F )T(k F )) FS /2|d(fe F )r r n 



7r|fcF±|£o-T(£ = kBT,k F ,a M ) ^n/ F s 

^ p& ^7 ^ ' F. T' <3) 



Vcos 2 9 k + sin 2 (0 fc - a M ) sin 2 9 k r ( £ = k ^ T -> fc F, «m) 
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where fepj_ is the component of k F projected onto the plane perpendicular to the field, and 
\k F \_\ = k F ^cos 2 9k + sin 2 (i^>fc — om) sin 2 Ok with k F = fe(cos(^fesin#fca+sin0fcsin0fcb+cos0jfcc) 
in polar coordinates [6]. a, b, c are orthogonal unit vectors fixed to the crystal axes. Here, 
we have assumed an isotropic spherical FS. The angle «m indicates the direction of the 
applied magnetic field rotated in a plane perpendicular to the c axis [6], and it is measured 
from the 0/t = direction (d-axis direction). r n is the scattering rate in the normal state. 
The coherence length is defined as £o = , v f ±/ttAq. We use the weak-coupling BCS ratio 
of Ao//cbT c = 2/1.13. The brackets (• • • )fs mean the FS integral with respect to k F like, 
(• • • }fs = (l/^o ) f dS F ■ ■ ■ /\v F (k F )\ = J dcftk J dOk sin^ • • • /(47r) with dS F being a FS area 
element and uq = f dS F /\v F (k F )\ being the density of states on the FS. The expression for 
r(e,k F ,a.M) is given in ref. [Bj. 

3. Results 

@We consider two Cooper pairing models on a spherical FS. One is the line-node s-wave pair: 
d(k F ) = \ cos (2 cf>k)\ sin 2 6k- The other is the d x 2_ y 2-w&ve pair: d(k F ) = cos(2(j)k) sin 2 9k- Both 



models have the same line nodes in the gap from the north pole of the Fermi sphere to the south 
one in the directions of fa = (1 + 2n)-7r/4 with n = 0, 1, 2, 3. On the other hand, the anti-node 
directions correspond to fa = rni/A. 

In Fig. [U we show the field-angle («m) dependence of pf , where each plot is normalized by 
its minimum pf m ; n and therefore the results do not depend on the parameter k-p^o or -fn/Ao- 
As seen in Fig. HJ^a), in the case of the line-node s-wave pair, pf exhibits its broad maximum 
when the magnetic field H is applied parallel to the gap-node direction (ajj ~ 0.8). It is noticed 
that pi has little temperature dependence. On the other hand, in the ci-wave case [Fig. QJb)], 
the sharp maximum appears when H is applied to the gap-node direction. The ratio of the 
maximum value to the minimum one is much larger than that of the s-wave case. In addition, 
the ratio increases with increasing the temperature T. 

The field-angle dependence of pf in the case of the constant r is displayed in Fig. [5J This 
corresponds to the case where we consider only the contribution of wo(^f) to pf neglecting the 
k-p dependence of the QP scattering rate r. The maximum of Pf («m) appears when the field is 
applied to the node direction, independent of whether the pair potential is the line-node s-wave 
or the d-wave. Such maximum direction «m in this case corresponds to the direction fa in the k 
space where uJo(kp) is small. The graph in Fig. [2] becomes flat if we consider an isotropic pairing 
d(kp) = const. 

Combined with our previous results for r(aM.) [S], the behavior of pf («m) can be understood 
as follows. In the line-node s-wave case, when the magnetic field is applied to the node direction, 
.T(ckm) takes the minimum [6] while Pf (om) due only to uJo(kp) exhibits its maximum [Fig. [2]. As 
a result, these two contributions to pf cancel each other and only the small and broad maximum 
appears [Fig. UK a)]. O n the other hand, in the d-wave case, r(ays) and pt(au) due to uo(kp) 
have their maximum in the gap-node direction in common. Therefore, the maximum of the 
field-angle dependence of pf stands out and the sharp maximum appears [Fig. [1Kb)]. 
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Figure 1. The field-angle («m) dependence of the flux-flow resistivity pf in the case of (a) the 
line-node s-wave pair and (b) the d-wave pair. The data are plotted for each temperature T. T c 
is the superconducting critical temperature. The vertical axis is normalized by minimum values 
Pf min for each temperature. 



4. Conclusion 

@We theoretically investigated the magnetic field-angle dependence of the flux-flow resistivity 
for the line-node s-wave pair and the d-wave pair on a spherical Fermi surface. We employed 
the quasiclassical approach and took into account a counterpart to the interlevel spacing of 
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Figure 2. The field-angle («m) dependence of the flux-flow resistivity pi in the case of r = const. 



the vortex bound states wo(fcF) as a contribution to the flux-flow resistivity, in addition to the 
quasiparticle scattering rate r(kp). The results show that the field-angle dependence of the 
flux-flow resistivity exhibits different behavior between the sign conserved line-node s-wave pair 
and the sign reversed d-wave one, irrespective of the fact that those pair potentials have the 
same amplitude anisotropy. Therefore, the field-angle dependence of the flux-flow resistivity can 
be a phase-sensitive probe if the flux-flow resistivity is observed, for example, by measuring the 
microwave surface impedance with changing the field-angle direction. 
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